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Abstract 



A Monte Carlo event generator is presented. An original algorithm is developed to 
simulate electron-positron scattering at energies and momentum transferred much more 
| than the electron mass. The first-order electroweak radiative corrections are included 

completely. Higher order corrections are taken into account by means of electron struc- 
ture functions. 



PACS: 12.20.-m Quantum electrodynamics, 12.20.Ds Specific calculations 
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^ ■ 1 Introduction 

a: 

The process of electron-positron (Bhabha) scattering was studied both theoretically and ex- 
perimentally since many years It has almost pure electrodynamical nature and therefore 
could be described with very high precision by means of perturbative QED. The process is 
commonly used at e + e~ colliders for luminosity measurements, because it has a large cross 
section and can be measured very accurately. The modern experimental technique of luminos- 
ity measurements reaches the one per mille level of accuracy, or even better, as at LEP1 ||. 
This is a challenge for the theory. At the Born and one-loop levels the process was investi- 
gated in detail (see papers H §, |5| and references therein), taking into account both QED and 
electroweak effects. The radiative corrections in the first order in the fine structure constant a 
become insufficient now, one has to take into account higher order effects. In order to meet the 
requirements of experiments one has also to implement the results of analytical calculations 
into a Monte Carlo event generator. 

In this note we present an event generatorr, based on the approach |]], how to merge the 
complete 0(a) result with the leading logarithmic corrections in higher orders. The structure 
of the paper is as follows. In the next section we give the master formula for description of 
large-angle Bhabha scattering and decompose it into 13 parts of different kinematics. The 
consequent structure of the Monte Carlo code is described in Sec. 3. The main options and 
parameters of the program are presented in Sec. 4. Numerical results and the precision achieved 
are discussed in the Conclusions. 



*on leave of absence from Joint Institute for Nuclear Research, Dubna, Russia. 
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2 The master formula 



The reaction 



e (pj) + e + (p 2 ) -> e (pi) + e + (j/ 2 ) + (717) 



will be considered in the centre-of-mass reference frame of the incoming particles. 
Let us start with the master formula in the form as given in paper |7j]: 
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Step functions represent any possible cuts on the phase space of the corresponding variables. 
By Ksv we denoted the so-called if-factor[], which comes from virtual and soft radiative 
corrections, and therefore it can be factorized at the Born cross section doo, 
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1 The last term in the square brackets of the expression for the if-factor in Ref. is incorrect. 
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The shifted (boosted) Born cross section (with vacuum polarization effect taken into account) 
reads 



d<r(zi,z 2 ) 
dipdc 



t 
a 



4a 2 
sa 2 

Re 



a 2 + zlil + c) 2 



;i -n(t)) 2 2z 2 (i -c) 2 

7 2 

-2 



+ 



z 2 (l-c) 2 + z 2 {l + c) 
|l-nm| 2 2a 2 



z 2 (l + c) 2 



[i -n(f»(i -n(s))*azi(i -c) 



-^sziYi(l - c), 

Zl + ^2 - (#L - 2k)c. 



2^12:2 



(4) 



For detailed notation and the derivation look in Ref. [/J. The current version of the code 
includes also the contributions of Z-exchange and Z-j interference as well as the relevant set 
of the first order weak radiative corrections to be described elsewhere. 

We quoted the complete expression; now we are going to discuss it and decompose into a 
form, suitable for an event generator. 

The first term of Eq. (0) is written in the form of the Drell-Yan cross section. The 
leading logarithmic corrections are accounted by means of the ©-functions, which are the 
kernel functions of the Dokshitzer-Gribov-Altarelli-Parisi-Lipatov evolution equations. The 
rest supplements the sub-leading terms, which come from the straightforward calculations in 
the 0(a) order. 

The non-singlet electron structure function is expanded in the series in a: 
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Here z means the energy fraction of an electron just before it emitted a collinear photon. The 
leading lograthmic corrections due to electron-positron pair production can be easily added 
within the same formalism |7|. H. 



2.1 The kinematical regions 

Our idea is to decompose the formula according to different types of the final state kinematics. 
Each part of the decomposition is a particular contribution to the total Bhabha cross section, 
and it can be measured, in principle, independently. So, as a result we have a sum of positive 
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quantities, which is important for Monte Carlo simulations. The 13 contributions are given 
below. 

1. The (quasi — ) elastic kinematics. 

Here we take into account the Born cross section with virtual loop corrections and soft 
photons. The energy of soft photons does not exceed Ae. The parameter A <C 1 is 
auxiliary, the final result (the total sum) should not depend on its value. That provides 
an additional check of our calculation. 



dipdc dipdc I 7r 7r 2 J 

2. k || p x : one (two) photons along the initial electron momentum. 

In this case we observe emission of a hard collinear photon inside a narrow cone along 
the direction of motion of the initial electron. The auxiliary parameter m e /e <C <^ 1 
defines the cone: kp 1 < 9 . We suppose that the parameter is less than the angular 
resolution of the detector. And therefore we are not going distinguish the situations of 
one or two photon emission in this cone. We just sum up the energies and momenta, if 
they are two. 
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The lowest value of z x to be defined from the conditions of particle registration. If y th is 
the threshold energy for electron registration, then 
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3. k || p 2 : one (two) photons along the initial positron momentum. 

This case is completely analogous to the previous one. 
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4. k || p[: one (two) photons along the final electron momentum. 
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The compensator for the final state radiation is 
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5. k || p 2 : one (two) photons along the final positron momentum. 
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6. fc a || p 1 , k b || p 2 : one photon along the initial electron momentum and one 
photon along the initial positron momentum. 
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7. fc a || p^, || p 2 : one photon along the final electron momentum and one 
photon along the final positron momentum. 
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9. k a || p 1 , k b || p' 2 : one photon along the initial electron momentum and one 
photon along the final positron momentum. 
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11. fc a || p 2 , fc b || p' 2 : one photon along the initial positron momentum and one 
photon along the final positron momentum. 
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12. || p 2 : both the final electron and positron go together back— to— back to a 
hard photon. 

This contribution works only if no any cut-off on acollinearity is imposed. Its kinematics 
is just the one of the process of annihilation into two photons (one of which is converted 
then into an electron-positron pair with low invariant mass). It gives a large logarithm, 
when the angle between the pair components is small. 
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where z denotes the energy fraction of the final electron, the new auxiliary parameter 
Ai bounds the energy of the hard photon for this contribution: 1 > uo > e(l — A x ). The 
opposite condition is to be implied in the last (WT) contribution in order to cancel out 
Ai. 
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13. One hard photon at large angle in respect to all other particle momenta 
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Letter denotes the general restrictions on the phase space: polar angles of the photon 
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cut-off on acollinearity for the final charged particles, if it is imposed. 
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The current version of the program contains the option to take into account the vacuum 
polarization affect in this contribution M. 

The complete final state kinematics can be defined in each contribution according to the 
general formulae given in Ref. J7[. 

3 Event generator structure 

An original algorithm for event generation was applied in the code. The steps are as follows. 

Step 1. At first we perform the numerical integration of the 13 contributions over the 
phase space, where only the most general cuts are applied: 

ai = IdT dTh * = 1 >---> 13 - ( 24 ) 

In this way we obtain the relative weights of the 13 contributions and at the same 
moment the absolute value of the total cross section: 

13 

Ui = —, o- tot = ^(Ti. (25) 

Ctot i= i 
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Step 2. The total ordered number of events to be generated N tot is shared between the 
13 contributions according their relative contributions to the total cross section. That 
is done as follows. The mail subroutine of the generator is called N tot times. Each time 
we choose one of the 13 kinematical regions according to the value of a random number 
r, by comparing it with the relative weights wf. 

j 3+1 

if ££4<r<]Tf4, then i = r, (26) 

k=l k=l 

if r < U\, then % = 1. 

Step 3. Now a set of kinematical variables v n for an event of the chosen contribution is 
generated. For each particular differential distribution dai/dTi we use a specific change 
of variables in order to make the distribution more flat. The weight of the event is 
defined by the formula: 

dai(v n ) N tot 
0.1 i cr tot 

Calculating the value of dcr, we apply the same cuts and conditions as while the numerical 
integration. 

Step 4. Because the differential distributions are rather complicated, we were not able 
to find a substitution to make them completely flat. So, the weights, obtained in the 
third step, can be different from unit. Here we do the following trick: generate a random 
number r and define the number of corresponding unweighted events as the integer part 
of w n + 1 — r. If the obtained number m n is more than 1, we use the rotation symmetry 
and distribute the events uniformly in the polar angle of the scattered electron. 

Step 5. At this step we can analyse the events generated in the previous step and apply 
additional cuts, if required. We can also record the events for future processing. 

Step 6. After we executed the steps 2-4 A^ot times, we can compare the results of 
numerical and Monte Carlo procedures. That provides a control of the technical precision 
of the code. Namely, we compare the total number of unweighted events with the ordered 
one: 

JVtot 

]T m n « N tot . (28) 

n=l 

At the same moment this means that the value of the cross section for the generated 
events is close to the one obtained by the numerical integration. Such a comparison 
is also done for each of 13 contributions separately. The technical precision can be 
improved by increasing of the total number of events and also by tuning parameters of 
the program for a concrete task. 



4 Flags and parameters 

The reading of flags and parameters is performed by means of the standard FFREAD subroutine, 
which is called from the PACKLIB JTOf. 
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The code contains several flags which can switch between different options in the physical 
base and in the generation procedure. Below we describe the most important flags. 



ICOR 
IC0R= 
IC0R= 
IC0R= 
IC0R= 

IORD 
I0RD= 
I0RD= 
I0RD= 

IVPDL 

IVP0L=0: 

IVP0L=1: 



=0: 
=1: 
=2: 
=3: 

=0: 
=1: 
=2: 



calculations only at the Born level 

Born + LLA corrections 

Born + LLA corrections + if-factor 

Born + LLA corrections + if-factor + large-angle photon 

calculations only at the Born level 

0(a) RC are taken into account 

higher order LLA corrections are included 

otQED is not running 

vacuum polarization by leptons and hadrons is accounted 



The $-meson contribution to the photon virtual propagator is realized in the program as a 
part of the vacuum polarization function. This option can be switched on/off by the flag 
IPHI=l/0. 

The main parameters to be set by user are: 
EB: the beam energy in GeV; 

TETN: the minimal electron scattering angle in radian; 
TETX: the maximal electron scattering angle in radian; 
TEPN: the minimal positron scattering angle in radian; 
TEPX: the maximal positron scattering angle in radian; 
NEVE: the number of events to be generated; 

TACO: the minimal allowed angle between the outgoing electron and positron in radian; 



In order to take into account vacuum polarization in the 13th contribution one has to 
set IEWT=1 and IVWT=1. This corrections is a part of the second order next-to-leading 
contributions (~ 0(a 2 L)), but it might be really important in the region close to resonance 
peaks. 



5 Conclusions 

The presented formulae are valid for the electron-positron colliders of moderately high energies 
below 3 GeV. In order to expand them for higher energies we take into account Z-boson 
exchange and the relevant electroweak radiative corrections |J, but the corresponding option 
of the code will be described elsewhere. 

In Table 1 we give the results for a rather simple configuration: beam energy 0.5 GeV, 
scattering angles for the both positron and electron lie between 15° and 165°, minimal angle 
between the final particles 30°, threshold for electron registration 50 MeV. For the Table we 
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Table 1: Results of numerical ca 


culations 




^Born 


°LLA 


°LLA+K 








without vacuum polarization 


Numer. int., mbarn 


13.501 


13.360 


13.023 


13.325 


13.303 


Monte Carlo, mbarn 


13.501 


13.359 


13.024 


13.324 


13.298 






-1.05 


-3.55 


-1.31 


-1.47 




with vacuum polarization 


Numer. int., mbarn 


13.749 


13.604 


13.261 


13.631 


13.638 


Monte Carlo, mbarn 


13.749 


13.603 


13.259 


13.627 


13.632 


5rc, % 


1.83 


0.76 


-1.78 


0.96 


1.01 



generated 10 8 unweighted events. The last line of the Table represents the relative difference 
of the cross section in the corresponding approximation in respect to the Born one: 

S RC = °i - aBorn 100%. (29) 

<^Born 

Further, in the Table we denote: cr Born is the pure Born level QED cross section; is the 
cross section with the O(a) LLA corrections; in 0lla+k ^ ne ®( a ) i^-factor due to soft and 
virtual photons is included; quantity cr^ includes the complete 0(a) set of corrections; the 
complete set of the first order corrections plus the LLA second order effects are presented 
in the last column. 

The resulting precision of the code for the description of large-angle Bhabha scattering in 
the typical conditions of electron-positron colliders of energy about a few GeVs is estimated to 
be 0.2% 0. The uncertainty can be decreased. In particular, an extended program for calcula- 



tions of second order next-to-leading RC to large-angle Bhabha scattering is in progress [|TT 
From the other hand, the uncertainty should be re-estimated taking into account the concrete 
experimental conditions and comparisons with other codes. 

Comparisons with other available Monte Carlo event generators for Bhabha scattering are 
in progress. The next step of the code development will be the explicit generation of radiated 
pairs. The third order leading logarithmic corrections will be also taken into account. The 
current version of the program does not include the strange effect of a diviation from the leading 
logarithmic approximation in the 0(a 2 L 2 ) order, which was recently found in Ref. [12 . 

The presented version of the code is dealing with the large-angle Bhabha scattering. That 
does not mean that the code can not be used for small-angles, but that it does not include 
the complete set of the second-order next-to-leading corrections, which are known only for 



the small-angle limiting case [O]. An extended version of the code, which will provide also 



the small-angle Bhabha event generation, to be described elsewhere. 
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